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Context. The weak gravitational lensing distortion of distant galaxy images ("sources") probes the projected large-scale matter distri- 
bution in the Universe. The availability of redshift information in galaxy surveys also allows us to recover the radial matter distribution 
to a certain degree. 

Aims. To improve the S/N in the mass mapping, we combine the lensing information with the spatial clustering of a population of 
/^^ ' galaxies that trace the matter density with a known galaxy bias. 

p^ ' Methods. We construct a minimum-variance estimator for the 3D matter density that incorporates the angular distribution of galaxy 

v,^ \ tracers, which are coarsely binned in redshift. Merely the second-order biasing of the tracers has to be known, which can in principle 

be self-consistently constrained in the data by lensing techniques. To study the new estimator, we generate a mock survey with galax- 
ies that trace the matter density with a Gaussian linear stochastic bias. 

Results. The filter smoothes and linearly mixes the individual lensing mass and tracer number density maps into a combined mass 
Q I map. The weighting in the mixing depends on the S/N of the individual maps and the correlation, R, of the matter and galaxy density. 

J-H , We down-weigh the influence of the dominant tracer number density by rescaling the noise covariance in the filter Even for a mod- 

f/j ■ erate mixing, the S/N in the mass map improves by a factor ~ 2, most strongly for z ^ 0.5. Moreover, the systematic offset between 

d ' a true and apparent mass peak distance ("z-shift bias") in a lensing-only map is eliminated, even for weak correlations of R ~ 0.4. 

The mixing, however, introduces a new noise component because of the stochasticity in the matter-tracer density relation. We give a 
prescription of the noise level and the average radial smoothing in the Gaussian regime. 

Conclusions. If the second-order bias of tracer galaxies can be determined from a lensing survey, the synergy technique potentially 
provides an option to improve the quality of a lensing 3D mass map. The detailed estimator bias and the noise, however, depends on 
Y-l ' the galaxy biasing mechanism, such that our Gaussian model may be too inaccurate at small angular scales. 

^SJ ' Key words. Gravitational lensing:weak - (Cosmology:) large-scale structure - (Cosmology:) dark matter - Methods: data analysis 

f^ 1. Introduction eral theory of relativity. Correlations between shear distortions 

CN of sources at different redshifts (cosmic shear tomography; Hu 

^~~l ' , . . , , . r- . „ ,,.,,, 2002) additionally constitute a promising probe for invesfigaf- 

• • The weak gravitational lensmg errect IS a well-established fool to ■ ^u i » j ■ e\u n • . j ue^ u i 

••^ • . r , ■ , , , ,, ■ ., ins the accelerated expansion or the Universe at redshirts below 

"^ mrer the protected, and to some degree also deproiected, large- i /c u \.i. i . i onin\ 

,..,., rr , i«r.>, n , z ~ 1 (Schrabbacket al. 2010). 

scale matter distabution (e.g. Tyson et al. 1990; Bonnet et al. 

»2] 1994; Small et al. 1995; Tyson & Fischer 1995; Squires et al. The coherent angular shear distortion pattern can be trans- 
it . 1996a,b; Fischer & Tyson 1997; Clowe et al. 1998; Abdelsalam lated into a map of projected matter fluctuations. Early non- 
et al. 1998; Fischer 1999; Clowe etal. 2000; Umetsu&Futamase parametric mapping algorithms, which were refined later on 
2000; Lombard! et al. 2000; King et al. 2002; Miyazaki et al. to obtain optimised methods for finite fields, achieved this 
2002; Dahle et al. 2003; Markevitch et al. 2004; Taylor et al. on the basis of a catalogue of source angular positions and 
2004; Jee et al. 2005; Dietrich et al. 2005; Clowe et al. 2006; ellipticities only (Kaiser & Squires 1993; Broadhurst et al. 
Bradac et al. 2006; Jee et al. 2007; Massey et al. 2007; Heymans 1995; Kaiser 1995; Schneider & Seitz 1995; Seitz & Schneider 
et al. 2008; Bradac et al. 2008; Holhjem et al. 2009; Gavazzi 1995; Schneider 1995; Seitz & Schneider 1996; Bartelmann 
et al. 2009; Okabe et al. 2010; Israel et al. 2010; Schirmer et al. et al. 1996; Squires & Kaiser 1996; Schneider 1996; Seitz & 
2011; Simon et al. 2012). The therein exploited coherent shear Schneider 1997; Seitz et al. 1998; Bridle et al. 1998; Seitz & 
distortions of distant galaxy images ("sources") are the result Schneider 2001; Saini & Raychaudhury 2001; Marshall et al. 
of the continuous deflection of light bundles by the intervening 2002; Starck et al. 2006; Marshall 2006; Deb et al. 2010; Deriaz 
fluctuations in the large-scale gravitational field (e.g. Schneider et al. 2012). With the advent of distance indicators of galax- 
2006a,b), most prominently detectable around galaxy clusters, ies in wide field galaxy surveys, the purely geometric relation 
The lensing distortions probe the total matter content in the between shear magnitude and source (and lens) distance is in- 
Universe, which makes them an excellent tool for studying the corporated into a new three-dimensional (3D) lensing algorithm 
dark matter component, an essential ingredient of the standard to also recover information on the radial distribution of matter 
cosmological ACDM model (e.g. Dodelson 2003), or for test- (Hu & Keeton 2002; Bacon & Taylor 2003; Simon et al. 2009; 
ing the dark matter hypothesis (e.g. Guzzo et al. 2008; Reyes VanderPIas et al. 2011; Leonard et al. 2012). So far, the best 
et al. 2010), which relies too a large extend on Einstein's gen- studied methodologies utilise linear inversion techniques, such 
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as Wiener filtering or a radial matter density eigenmode decom- 
position with a suppression of low signal-to-noise (S/N) modes. 
Owing to the relatively sparse and noisy sampling of the sur- 
vey area with background sources, however, the resulting maps 
are usually very noisy, and significant detections are basically 
restricted to mass peaks of a galaxy cluster scale with only mod- 
erate redshift accuracy; the mass detection limit and redshift res- 
olution declines quickly with peak distance (Simon et al. 2009, 
2012). Moreover, because of the low S/N in mapping applica- 
tions, the linear inversion applies a radial smoothing with a broad 
smoothing kernel that (a) smears out localised peaks in radial di- 
rection and (b) biases the peak distances ("z-shift bias"; Simon 
et al. 2009), which potentially renders the resulting maps hard 
to interpret. To attain more realistic 3D maps, the radial elon- 
gation of peaks inside the map can be mended by regularising 
the inversion (Leonard et al. 2012), or by finding the maximum- 
likelihood (ML) positions of one or a few individual mass peaks 
along the line-of-sight (l.o.s.) given the radial smoothing kernel 
and radial density profile in the map (Simon et al. 2012). But 
this does not alleviate the principle problem of noisy maps and 
inaccurate peak distances. It merely provides more realistic ML 
estimators for the 3D mass map. Moreover, the noise properties 
of the maps are likely to be complex in regularised, non-linear 
methods. 

On the other hand, galaxy positions themselves are also trac- 
ers of the 3D matter density field and could therefore be em- 
ployed to add extra information to the matter density maps that 
are obtained from 3D lensing. There are two complications here, 
however: (i) galaxies trace the matter density field only up to a 
systematic mismatch, which is generally dubbed "galaxy bias" 
(e.g. Weinberg et al. 2004), and (ii) a sampling by galaxy posi- 
tions is affected by shot-noise (e.g. Martinez & Saar 2002). 

Pen (2004) discussed the cross-correlation of galaxy posi- 
tions, endowed with highly accurate spectroscopic redshifts, and 
two-dimensional (2D) lensing matter maps to improve the con- 
straints on the matter power spectrum and facilitate the depro- 
jection of lensing maps. The strategy of this paper is to refine 
the minimum-variance estimator in Simon et al. (2009) for the 
3D matter density by adding the galaxy clustering information to 
the map making process. In contrast to Pen (2004), galaxies do 
not need to be equipped with accurate redshifts, however, broad 
redshift slices, as expected from photometric redshift surveys, 
suffice. Since the minimum-variance estimators (Zaroubi et al. 
1995) only require 2"''-order statistics of the input data to be 
specified, only the 2'"'-order bias parameters (Gaussian biasing 
or linear stochastic bias; Pen 1998) of the galaxy tracers have to 
be known. The galaxy bias as function of scale and redshift could 
in principle be acquired in a self-consistent approach from the 
data by using lensing techniques (Schneider 1998; van Waerbeke 
1998; Pen et al. 2003; Fan 2003; Julio et al. 2012; Simon 2012), 
or with lesser certainty from simulations (Weinberg et al. 2004), 
such that we henceforth assume that it is basically known. As to 
(ii), the galaxy sampling shot-noise is taken care of by the galaxy 
noise covariance within the minimum-variance estimator. 

Thus, we pursue here the idea to combine 3D lensing and 
galaxy clustering information to map the 3D matter distribu- 
tion. Towards this goal. Sect. 2 summarises the formalism of 
Simon et al. (2009) to obtain the matter density on a set of 
lens planes based upon a set of source samples binned in red- 
shift. Furthermore, this section lays out a new formalism to ob- 
tain the galaxy number density on the lens planes for a given 
set of galaxy tracers, binned in redshift slices, with known lin- 
ear stochastic bias. A synergy of both estimators is achieved by 



Sect. 3. The benefits of a synergy are explored in Sect. 4 using 
simulated data. Sect. 5 finishes with a summary and conclusions. 

2. Independent reconstructions 

2. 1 . Matter density on lens planes 

Here, we briefly summarise the formalism already presented in 
Simon et al. (2009). We adopt the exact notation that is employed 
therein. For more details, we refer the reader to this paper. 

We split the source catalogue into / = 1 . . .A^z subsamples, 
all of which with known redshift probability distribution (p.d.f.). 
The complex ellipticities (Bartelmann & Schneider 2001) of the 
sources belonging to the /th subsample are binned on a 2D grid 
that covers the field-of-view of the survey area. This ellipticity 
grid is denoted by the vector e*'', whose elements are the sorted 
pixel values of the grid. Every subsample uses the same grid 
geometry. 

The paper assumes that the weak lensing approximation is 
accurate enough for the lensing catalogue on the whole, i.e., the 
complex ellipticities, 6\ of individual sources are unbiased esti- 
mators of the shear, y, for the given source redshift in the l.o.s. 
direction 6; of the source. 



<£^> = (7 + e^) = y , 



(1) 



where e' denotes the intrinsic (unlensed) complex ellipticity of a 
source image. Moreover, we assume a flat-sky with a Cartesian 
coordinate frame. 

We slice the light-cone volume, for which the matter dis- 
tribution is reconstructed, into N\p slices. Within the slices, the 
matter density contrast is constant along the l.o.s. as approxima- 
tion. The solid angle about a l.o.s. direction corresponds to one 
pixel of the grid. Thus, the matter density field inside a slice is 
fully described by the angular distribution of mean density con- 
trasts on a plane (lens plane) and the width of the slice. The mat- 
ter density contrast on a lens plane, ff^, is binned with the same 
grid as the source ellipticities. We represent the grids, e''' and 
tfJn, as vectors of equally ordered pixel values. We refer to a par- 
ticular pixel by 6^(0/), where 9j is the direction of the pixel on 
the sky. Therefore, our algorithm represents the 3D-matter den- 
sity contrast as an approximation by a discrete set of lens planes, 
which numerically limits the radial resolution, and a discrete set 
of pixels on the sky, limiting the angular resolution. 

The complete sets of ellipticity planes and lens planes are 
combined inside vectors of grids. 






c(A'z) 



]. 



c(A'ip)l 



(2) 
(3) 



respectively. The brackets, grouping together the vector argu- 
ments, should be understood as big vectors that are obtained by 
piling up all embraced vectors on top of each other. 

In the weak lensing regime, the (pixelised) lensing conver- 
gence K^'XOj) is, in the lowest-order Born approximation, the 
weighed projection of the density contrast on the lens planes: 

" Mp Wip 

^= J]QiiS'-^,---,J]QNj6'^ =:Q.5™, (4) 



where the coefficients Qij express the response of the /th con- 
vergence plane *:''' to the density contrast in the y'th lens plane, 
namely 



Qu 



3Qm P 

2Dl 4 



dx 



a(x) 



(5) 
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where 

^W= I dx "^'^,^,C' pf(x'), (6) 

Px (x) denotes the p.d.f. of sources in the /th source subsam- 
ple, and L^j,;!^ j+i [ sets the comoving radial boundaries of the jth 
matter slice. By Dfj := c/Hq, we mean the Hubble radius and by 
/kC^) the (comoving) angular diameter distance. The projection 
from a grid vector in "<J,„-space" to a grid vector in "iif-space" is 
hence denoted by the operator Q that is acting on 6m- 

The next step connects the convergence planes k to the shear 
planes by a convolution of the lensing convergence on the grid. 



Y ~ ^yK^ 5 • ■ ■ ! r y/f^ — . VyiiK , 



(7) 



which introduces the operator Py^ (Hu & Keeton 2002) to map 
jic''' to the corresponding shear plane 7*''. In this sense, Py^ per- 
forms a linear transformation from k to 7-space. 

Using this compact notation, we express the linear relation 
between the matter density (contrast) on the lens planes and the 
observed, binned, ellipticity planes as: 



e = Py;,Qtf,n + riy 



(8) 



Here, an additional vector ny denotes the binned intrinsic ellip- 
ticties of the sources of all source subsamples. In the language 
of lensing, we consider this the noise term that dilutes the shear 
signal Py^OS^. 

For the scope of this paper, possible correlations between 
shear and intrinsic shapes are ignored (Hirata & Seljak 2004). A 
minimum-variance estimator of dm in Eq. (8) is then 



■'m,est 



S,Q'P;,(N;'Py,QS,Q'P;, + al) 



n: 



(9) 



As only input, the minimum-variance filter requires the signal 
covariance 8,5 - {dmS\^), which specifies the presumed two- 
point correlation between pixel values of S^iO) on the lens 
plane(s), and the noise covariance Ny - {nyuL), which quan- 
tifies the shear pixel noise variance and the correlation of noise 
between different pixels. Pixels that contain no sources have infi- 
nite noise. For the prior signal covariance, correlations between 
pixels that belong to different lens planes are set to zero. We note 
here that the signal covariance does not need to be the true sig- 
nal covariance in the data, although the reconstruction may be 
suboptimal when it is not. 

The signal covariance determines the degree of smoothing 
in the 3D map. The smoothing is uniquely defined by the linear 
transformation 



B, 



5 ■- 



(ff1 + S,N,')"' S,5N,' ; N,' := Q'Pj.N^'Py.Q , 



(10) 



and can as such be utilised for a comparison of the map (5,n,est 
to a theoretical matter distribution dmxh by Ba(5m,th (Simon et al. 
2012). The radial smoothing is characterised by a radial point- 
spread function (p.s.f.) of the filter. After smoothing with the 
radial p.s.f., a peak in the true matter distribution (5m,th does 
on average not necessarily peak at the same distance as in the 
smoothed map, which gives rise to the so-called redshift bias or 
"z-bias". Inside the filter, the constant a e [0, 1] tunes the level 
of smoothing by rescaling the noise covariance. 

From a practical point of view (Appendix B of Simon et al. 
2009), the Wiener filter consists of a series of linear operators 
that is applied step-by-step from the right to left to the grids. 
Within this process, the signal covariance Sg is a convolution or. 



equivalently, a multiplication in Fourier space of Fourier modes 
f(() of the /th lens plane with the angular signal power spectrum 
P^{(), which is implicitly defined by 



d'ms'rlit')) = i2nf5^{t + t')P7m) 



m. 



(11) 



We approximate the power spectrum by using Limber's equation 
(Kaiser 1992) in Fourier space. 









(x)r''\fK(xy'^l ' 



(12) 



where Axi := Xi+i ~ Xi^ ^(0 is the Fourier transform of the pixel 
window function, and P3d{k,x) is the 3D matter-density power 
spectrum at radial distance x for wave-number k (e.g.. Smith 
et al. 2003), 6d{x) is Dirac's delta function. We denote the 
Fourier transforms of flat fields f{0) on the sky by /(/"), defined 
by 



f(0) 






-J' 



-W{ 



f{t)c^"'' ; fit) - d^0/(0)e-"'' . 



(13) 



2.2. Galaxy numbers densities on lens planes 

To improve the information in the 3D matter map and to pos- 
sibly alleviate the z-shift bias, we add the information gained 
from galaxy positions, which also probe the matter distribution 
("tracers"). 

In this section, however, we first visit the problem of map- 
ping the spatial galaxy number densities. For this purpose, we 
estimate the number density of galaxies projected onto the pre- 
viously defined lens planes. Hence, we slice the full true 3D 
galaxy distribution into A^ip distance slices with distance limits 
[A'/'A'!+i[- The galaxies are counted within each slice and angu- 
lar grid pixel of solid angle A„. Thereby, we receive the galaxy 
number density tig (6j) - N'^'^{0j)IA^ in l.o.s. direction 6j of the 
/th slice, where N^'^Oj) is the number of counted galaxies. We 
compile the galaxy number density values inside a grid vector 
M„ , and we then arrange all grids inside a vector of grids 



n,^\£\. 



„Wp)l 



(14) 



This is the number density distribution of galaxies that the 
following scheme seeks to recover from a galaxy sample with 
inaccurate distance information. Towards this goal, we split the 
observed galaxy sample as to their redshift estimators, Zest 6 
[z(Xi),z(Xi+i)L into A^ip subsamples with known radial p.d.f. 
Pf (x). By projecting the /th sample onto a 2D grid on the sky, 
one obtains the observed number density distribution 

A'lp Mp 

r^iOk) = 2 f'-^^^(^k)Pu nfiOtd =: ^ G,/0t)«^^\0t) , (15) 
where /mask 6 {0,1) flags mask pixels (= for mask), and 
p,,:= r" Axpfix) (16) 

is the probability to find a galaxy with;if 6 [;ify,A'y+i[ anywhere 
inside the grid rj"^ . Owing to the redshift errors, the observed 
distribution on the lens planes, ri„ (Ok), does not exactly match 
the true distribution n^\0k). Therefore, < Gij{0k) < 1 denotes 
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the fraction of galaxies on the jth lens plane that are, on average, 
mapped onto the grid t;! . Because of masking, the total number 

of galaxies is not necessarily conserved, i.e., 2,j'i Gij{0k) + 1. 
By a proper arrangement of the elements G\j{Ok) inside a matrix 
G, the effect of Gij{0k) on the entire 3D grid /ig can be written as 



% 



G«„ + ( 



where 



-Wi 



.4"'1 



(17) 



(18) 



We presume that galaxies sample an underlying smooth galaxy 
number density by a discrete Poisson process (e.g., Martinez & 
Saar 2002). Therefore, the observable galaxy counts sample the 
underlying galaxy number density Kg merely up to shot-noise, 
which is here formally expressed by the noise component ^g. 

Now, by analogy with the matter density d^, we can find 
an minimum-variance filter to estimate the true distribution of 
galaxies on the lens planes, namely 



*g,est 



= s,G' n:1gs,g' + 



.1)- 



n; 



%• 



(19) 



As before, Sg - («g«g) is the signal covariance, which is the an- 
gular clustering two-point correlation function of the galaxies on 
the lens planes, and Ng = (^g^L) denotes the shot-noise covari- 
ance. The degree of smoothing by the Wiener filter is tunable by 
using y6 € [0, 1], which does not need to equal the parameter a in 
Eq. (9). For the Poisson shot-noise covariance, we adopt a diag- 
onal noise covariance, [Ng],j = for / ^ j, with [Ng]„ = 77! '(61) 
for unmasked grid pixels 0i, and infinite noise otherwise. ' By 
fig {61), we denote the estimated mean number density of galax- 
ies in pixel 0i of the A:th subsample (see next section). 

Like for the matter density Wiener filter, a practical imple- 
mentation of the Wiener filter in Eq. (19) consists of a series 
of linear operations applied to //g. The effect of Sg is to multi- 
ply every angular mode f({) of the /th lens plane with the prior 
galaxy power spectrum Pg ({), which we define relative to the 
matter power spectrum using the galaxy bias factor /7*'*(^) > 
(e.g., Tegmark & Peebles 1998), 



(hfiOhfie')) = (2nf6t,{{ + r) [nfb''\\e\)f pfm 



(20) 



rt'w 



^(0 



where «g denotes the true mean number density of galaxies on 
the (th lens plane. For this definition of the bias factor, the shot- 
noise contribution to the galaxy power spectrum is excluded as 
it is already accounted for in 0g. 

The angular bias factor b^'\{) is related to the 3D bias factor 
b{k,x), where k is the comoving 3D wave-number, by a projec- 
tion that is approximated by Limber's equation. 



[b^'\{)fpfi{) 



(21) 






• X 



P^ (0 is given by Eq. (12). For this approximation, we assume 
that the number density of galaxies stays constant inside the 
slice. 



' The Wiener iilter in the given form requires the inverse noise covari- 
ance, such that elements with infinite noise on the diagonal are zero. 



2.3. True mean galaxy numbers 

The true galaxy number densities hg in Eq. (20) have to be de- 
rived from the data itself. For an estimator of «„ , we go back to 
Eq. (17), which relates the observed number of galaxies, i/g, to 
the true number on the lens planes, «g. For an statistical ensem- 
ble average of this relation, we expect 



^g := (% 



= G«c. 



(22) 



;;('■)/ 



-.(•'> . 



wherein, however, all elements fig (Oj) equal one number ng ow- 
ing to the statistical homogeneity of the galaxy number density 
fields, hence 



A'lp 

f]'i'i0k)^J]G,j{ek)fi 



U) 



(23) 



Summing over all pixels, A^p in total, of the /th tracer sample 
yields 



^-^Z€(«^) = Z^'^^. 



(» 



P A = l 



J=i 



where 



G(/ 



Wp 

JVp'^G,v(0*) 

A-=l 



(24) 



(25) 



averages G^ over the area of the grid. Inverting the former equa- 
tion, gives 

N,„ _ N„ 



-JO 



^Z° 



X. 



;(j) 



X 



J) 



1 



J]r^f{e,). 



(26) 



P t= 



/t=i 



;('■) 



For an unbiased estimator of ?(„ on the r.h.s., we insert the ob 



/')/ 



J'h 



served galaxy number densities, i.e., rjg (0^) = rjg (0^). The value 

of ffg {0k), which is utilised for the noise covariance Ng in the 
foregoing section, is computed from Eq. (23) and the estimated 
frf. 

In the simple case of negligible redshift errors, we find 
Giji0k) = Sfjmask{0k), whcrc Sf- dcnotcs the Kronecker symbol. 
In this case, we consequently find 

[G h-S^^^ (27) 

for A^ being the number of unmasked pixels. Moreover, we find 
fig^ - A^pM'V(A^O) for a number Ng'' of galaxies within the /th 
subsample and a survey area Q. Thus, the galaxy number density 



j(0 



Ng /Q. is scaled up by Np/N to account for the mask. 



,('■) 



However, the estimator in Eq. (26) has one caveat since rig 
are basically a convolution of hg with the redshift error of galax- 
ies. A deconvolution through G possibly results in oscillating 
and negative values for fig . We therefore regularise Eq. (26) by 
a constrained deconvolution fig that maximises the likelihood 

In X(ngl^g) = (Gng - ^g)' N^ ' (G«g - ^g) (28) 

under the condition that «g > for all /. We determine this so- 
lution numerically. The additional covariance N,, can be used to 
give different weights to the observed ^g values, for example by 
weighing with the number of galaxies in each galaxy sample in 
order to account for the galaxy shot-noise. For equal weights, we 



simply set N;, 



1. 
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3. Combined reconstruction 

3.1. Minimum-variance estimator 

Up to here, we have considered the galaxy number density and 
matter density fields separately. However, //g contains informa- 
tion about 6m and vice versa, as galaxies trace the matter distri- 
bution to a certain degree. On a statistical level, this relation is 
reflected by a non-vanishing cross-correlation 



Sfe = <(Sm«' > 



(29) 



for pairs of pixels on the same lens plane, which has not entered 
our formalism so far. Slices are thought to be wide enough, such 
that correlations between pixels belonging to different lens plane 
are negligible. 

We combine the 6^- and «g-grids inside one new vector 



S := [(Jm,«g] , 



(30) 



which relates to the observed shear and tracer number density 
grids. 



according to the Eqs. (8) and (17) via 



Py^Qdm, G«g 



+ n -: Hs + n . 



where the combined noise vector is 

n :- \ny,<pJ . 

In this compact notation, the action of a matrix 



A = 



All Ai2 



A21 A22 

on a product vector v - [vi , V2] is defined as 
Av :- AiiVi -t- A12V2, A21V1 + A22V2 . 
In this sense, the projection matrix R is 

Py.Q 



R = 







G 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



Following the usual assumptions of a minimum-variance fil- 
ter, the optimal filter for estimating s from d in this combined 
problem is 



which uses as short-hand notations 



S = 



S, S 



<sg 



^Sg ^<^ 



; N„« = 



aHy 
/3Ng 



R"' = 



Q'p;. 

G' 



(37) 



.(38) 



The galaxy shot-noise tp^ and the intrinsic ellipticities of the 
sources, which are comprised in «y, are assumed to be uncorre- 
cted. By choosing different tuning parameters a ^ /?, the impact 
of the Wiener smoothing can be adjusted independently for the 
matter and galaxy map. 

The novelty of the combined reconstruction is that tracer 
number and matter density maps will exchange information, if 
the cross-correlation matrix S^g is non-vanishing. In a practical 
implementation of the filter (37), we apply, as before, step by 
step linear operations to the grids stored inside d. As with the 



previous operators Ss and Sg, the application of S^g amounts to 
a multiplication of angular grid modes with the cross-correlation 
power spectrum P''^\{), determined by 

«®(f)5«(r)\ - (2nf6^{t + {') hfr^'\\t\)b^H\t\)Pf{\t\) . (39) 



=P5g(l) 



d(0 



We define P\ {£) with respect to the matter power spectrum 
Pg (£) by employing the galaxy-matter cross-correlation factor 
r'-'\{) (Tegmark & Peebles 1998). The angular function /''(O is 
approximately related to the 3D correlation factor r{k,x) accord- 
ing to 



(,('■) 



r^'>i{W'>{OP!i{) 



,(0/ 



(40) 



\m\- 



2 p.^ 



dx 



[^Xi]-Jx, [/kOt)]- 



■r{ki,x)b{ke,x)P3d(kf,x) , 



where k( := ^//kOt)- 

To understand the mode of operation of the minimum- 
variance filter in Eq. (37), it is instructive to recast it into the 
mathematically equivalent form 



*est = (1 + SN,„'j"' SN,„'^ X Ns„^R^N-' d , 

Step-l 



(41) 



Step-2 



where N. ' := R'N iR denotes the inverse noise co variance of 
the matter and tracer density maps after the application of the un- 
biased estimator in Step-l.^ Step-l involves no Wiener smooth- 
ing, hence matrix S, to construct the maps. As this is usually too 
noisy, we apply an additional smoothing to these maps by virtue 
of the Wiener filter in Step-2. This filter linearly combines, aver- 
ages pixel values in the maps based upon the expected S/N in the 
unbiased maps. It is this Step-2, the analogue of the matrix Bg in 
Sect. 2.1, that introduces biases into the maps, especially through 
a radial smoothing. Moreover, only Step-2 formally mixes pixels 
from the mass map and the tracer number density map by means 
of the off-diagonal matrix S^g. Therefore, Step-l makes inde- 
pendent mass and tracer maps that are only later on combined 
in Step-2 according to our prior knowledge of the correlation of 
both. Setting a - p -Q results in a unity matrix for Step-2, i.e., 
no smoothing. 



3.2. Fourier space representation 

For shear and galaxy number noise homogeneous over infinite 
grids with no gaps, the estimator in Eq. (37) takes a simple form 
in Fourier space. Under this idealistic conditions, the angular 
modes 



S(f) = [j„(f),«g(f)] 



(42) 



are only linear functions of the //g- and y-modes of the same I; 
there is no mixing between modes of different E. Therefore, a 
reconstruction is then done most easily in Fourier space by 



7est(^) = SfRj(RfSfR; + H^pY'dit) , 



(43) 



^ Our implementation of the estimator explicitly assumes that both 
the matter and the tracer density contrast vanish when averaged over 
each lens plane to lift the mass-sheet degeneracy. 
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where d(() - Wit),]] (()\ are the observable input grids. The 
covariance matrix of the (homogeneous) noise is 



N» 



'yS 






ao-; , 



,m 



(Mp) 



(44) 



where ris is the mean source number density of the /th source 
sample (out of in total A^,), crl = (e'[e']*) is their intrinsic shape 
noise variance, and fig is the Poisson shot-noise power (white 
noise). Possible noise contributions owing to intrinsic align- 
ments of sources are ignored here, hence N,,/? has no ofF-diagonal 
elements. Furthermore, one has^ 



5 lD(e)Q 
^'^[ G 



D*(e)Q' 

g' 



(45) 



where D(() = (/e* (Kaiser & Squires 1993); for ^ = 0, we set 
D{{) = 0. The signal covariance is 






S,g(^) 



with 
S,(^ - diag{p^'>(^). 



,Pf{0, 



r-Wp) 



(O) 



Sg(^ = diag{...,[«l,'V)(^)]V®(^),...} . 



(46) 

(47) 
(48) 
(49) 



Because of the diagonal structure of the last three matrices, 
the matrix Sf acting on a vector ?(/") actually only mixes the 
matter and tracer density modes from the same lens plane. Thus, 
rearranging the grids inside v(f ), and pairing together matter and 
tracer density modes of the same lens plane, renders St a diago- 
nal block matrix 



S, = diag{s<",...,Sr^'} 
with the 2 x 2-blocks 



c(0_( 1 nfr^H^b^He)]^p(i).f. 

^ I «®r»(^)M')(^) [«®^®(^)]2 1'''^''^^^ 



(50) 



(51) 



4. 1 . Fiducial survey parameters 

For the fiducial survey, we split the source galaxy catalogue into 
A^, = 20 equally sized redshift slices of width Az = 0.1, span- 
ning the redshift range z - . . . 2.0. For the sources, we neglect 
the effect of redshift errors greater than the width of the redshift 
slices, such that the true p.d.f. pi {z) of sources of the /th slice 
is well-described within z G [z,, z,+i] by the p.d.f. of redshift es- 
timates of the full sample. 



pr{z) oc z^ exp - — 



-i\ 



(52) 



where zo = 0.57, A - 1.5, andz; - (i- l)Az. We represent the re- 
construction volume by A^ip = 10 lens planes between z = and 
z = 2, centred within slices of moderate width Azip = 0.2. The 
total number density of sources on the sky is h - 30 arcmin"^ 
with an intrinsic shape noise of cr^ - 0.3. 

To support the matter density reconstruction, we include 
fiducial galaxy tracers with known galaxy bias. For simplicity, 
their p.d.f. of redshift estimates is identical to Pziz). Contrary 
to the sources, however, we now also emulate the effect of red- 
shift errors by adopting a root-mean-square (r.m.s.) accuracy 
of o-~{z) - 0.04(1 +z) (Gaussian errors), which is built into G 
in Eq. (15). The slicing scheme for the tracers is equivalent 
to the scheme of the sources. From this, we compute the av- 



erage number density of tracers «g and ?7g from Eq. (23) for 

each redshift slice, and the true redshift distributions Pf (z) by 
piecewise convolving the p.d.f. Pz(z) with a Gaussian kernel of 
r.m.s. <T-(z)- For low redshifts, we have ff^ a; «g , but we find 
differences at higher redshifts where Azip ~ cr.. For the fiducial 
survey, we reduce the total number of tracers to 10 percent of 
the sources, Wg = 3 arcmin"^, since a reconstruction realistically 
requires a specifically selected tracer population for an accu- 
rately known bias. Here, the tracers are clustered as the mat- 
ter, b{{) = 1 for all redshifts, but, more relevantly, we assume a 
slight stochasticity on all scales, namely r{{) = 0.8. A high cor- 
relation, r > 0.5, for various galaxy populations is expected from 
theoretical models (e.g., Guzik & Seljak 2001) and observed for 
some cases (Hoekstra et al. 2002; Simon et al. 2007; Julio et al. 
2012). As model for the cosmic average 3D matter power spec- 
trum P3d{k,x)7 we devise the prescription by Smith et al. (2003), 
which is used for the signal matrix S^. 



on the diagonal. This structure is useful when implementing the 
estimator in practise. Clearly, modes will not aff'ect each other 
when r^'\£) = with no improvement by a synergy of lensing 
and galaxy tracers. 



4. Mock survey study 

We move on to quantify the impact of a joint reconstruction 
by considering an idealised survey with homogeneous noise 
and a G that is independent of the pixel position. As a fidu- 
cial cosmology, this paper uses a standard flat ACDM model 
with matter density parameter Q,n - 0.27, of which baryons are 
D.h = 0.046, and a Hubble parameter Ho - h 100 km s"' Mpc"' 
with h — 0.704. The normalisation of the matter fluctuations 
within a sphere of radius 8 A^'Mpc at redshift zero is erg - 0.81. 
For the spectral index of the primordial matter power spectrum, 
we use «s = 0.96. 



4.2. Signal-to-noise improvement 

With the estimator (43) at hand, we forecast the S/N of the matter 
and tracer number density modes as function of angular wave- 
number i. To this end, we compare, by analogy to Simon et al. 
(2009), the cosmic average power spectrum of the reconstructed 
matter and galaxy number density modes on the lens planes. 



Ps(^) := W,SfW; , 

to the predicted noise in their reconstruction, 

PM := W.XfWj , 

where we use 

Xf := (RfN^^Rff' ; W, := Sf (Sf + Xf)"' . 



(53) 



(54) 



(55) 



Here, G does not depend on ft. 



In this setup, the Wiener filter Wf devises the actual signal power 
S[ present in the data. As pointed out earlier, this is not a neces- 
sity, but is required for an optimal minimum-variance filter. 
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Fig. 1. S/N of density modes for tfie matter density maps (left panel) and galaxy tracer density maps (right panel). The lines assume cosmic-average 
fluctuations in the density fields and b{£) = 1 for the tracers. DilFerent lines styles correspond to different lens plane redshifts: z = 0.1 (solid), 
z = 0.3 (dashed), z = 0.5 (dashed-dotted), z = 0.7 (dotted), z = 0.9 (dashed-dotted-dotted). The black thin lines do not employ any smoothing, 
while red thins lines use smoothing; both use no mixing of lensing and galaxy clustering information. In the right panel, red and black thin lines 
coincide. Thick lines depict values in a Wiener smoothed map (a = 0.01, yS = 0.1) and a mixing with r(C) = 0.8 for all redshifts. 



For a - /3 - (no smoothing), the noise covariance is 
Pn(^) - X(. For the following, we set a = 0.01 to be consistent 
with Fig. 1 1 of Simon et al. (2009) for a lensing-only recon- 
struction with a transverse filter The parameter a must not be 
too close to unity, as this results in a too much radial smooth- 
ing, which moves basically all mass peaks to the middle of the 
reconstruction volume (no radial information). In order to give 
the tracers not too much weight in the joint reconstruction, we 
set yS = 0. 1 . A parameter /3 smaller than a is a means to down- 
weigh the impact of the tracers in the joint reconstruction. By our 
choice of parameters, the Wiener filter assumes that the galaxy 
shot-noise is larger by a factor of '\JpJa a; 3.2 with respect to 
the shape noise, thereby giving less weight to the galaxy clus- 
tering information. In addition, a = 0.01 utilises a Wiener filter 
that assumes shape noise smaller by a factor of ^/a = 0.1 than in 
the data and thus applies less smoothing to the maps. Our goal 
is to aid the 3D mass mapping a little by adding some galaxy 
tracer information, rather than being completely dominated by 
the tracers, which usually probe the matter density with consid- 
erably larger S/N compared to lensing - at least, if the galaxy 
bias is known; compare the thin black lines in the left to the red 
lines in the right panel of Fig. 1 . 

Fig. 1 depicts the S/N of the lens plane density modes as 
function of angular scale and lens plane redshift. In the left 
panel, we have the matter density modes, in the right panel the 
tracer number density modes. Different line styles correspond to 
different lens planes, the thin lines to reconstructions with mode 
mixing switched off, i.e., r{() = 0, and the thick lines to the 
joint reconstruction. Clearly, as discussed in length in the recent 
literature (Hu & Keeton 2002; Simon et al. 2009; VanderPlas 
et al. 201 1), a lensing-only map absolutely requires some radial 
smoothing, which is seen here by comparing the low S/N of the 
thin black lines to the boosted S/N in the red lines. Adding mode 
mixing with r{{) = 0.8 and the previous tuning parameters, fur- 



ther boosts the S/N in the mass map; here by a moderate factor 
of about two (red lines vs. thick black lines). An exception are 
angular scales of roughly one arcmin or smaller, where the shot- 
noise of the tracers starts to dominate their clustering signal. 

The impact of mixing on the S/N of the tracer number den- 
sity maps (right panel) is small, most prominently on the small 
angular scales. This changes slightly, if we choose an even larger 
tuning parameter /3 (not shown): a larger /3 scales up the shot- 
noise matrix of the tracers inside the Wiener filter, attributing 
even more weight to the lensing data in the joint reconstruction. 
As the S/N of the tracers in the data is actually higher than that of 
the shear, this will result in a decreased S/N for the galaxy num- 
ber density maps in comparison to a reconstruction with no mix- 
ing; the joint reconstruction is not optimal as to the map noise. 

We also observe here that for a cosmic average power in the 
density fields, the reconstruction algorithm recovers the density 
modes only with a poor S/N < 1 on sub-degree scales for the 
fiducial survey. This reiterates the past observations that map- 
ping with lensing, or in a combined reconstruction still domi- 
nated by lensing, may only be feasible for mass peaks largely 
exceeding the cosmic average power spectrum, such as galaxy 
clusters. 



4.3. Radial point spread function 

We now explore the impact on the radial p.s.f. in the matter maps 
by mixing 3D lensing information with the number densities of 
tracers. The radial p.s.f. is the average sight-line profile in the 
smoothed matter map of a single mass peak that is situated at 
a lens plane in the un-smoothed matter field. Ideally, the p.s.f 
spikes at the true mass peak redshift (no z-shift bias). We as- 
sume a homogeneous survey, such that the choice of the l.o.s. 
direction 6 is irrelevant. We hence arbitrarily pick = as ref- 
erence direction and omit the pixel index in the following. 
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Fig. 2. Tlie radial p.s.f. in ttie smootiied mass maps for tiie case r{C) = 0.8 (right panel) and no mixing of lensing and galaxy clustering {r{{) = 0; 
left panel). The details of the fiducial survey are found in Sect. 4.1. Both panels adopt a = 0.01 andjS = 0.1. The redshifts of the density peaks in 
the un-smoothed maps are the small number labels, which are only shown up to z = 1.3. The map pixel size is ©s = 1 arcmin. 



We consider a singular test density peak with profile 
•JmC^) = Ap5u(x ~Xi) in the un-smoothed map. It is located at 
the distance^/ :- (xi + Xi+i)l'^ of the /th lens plane. For circular 
pixels with angular radius ©s, the pixel value of this peak is in 
the un-smoothed map 



F(A6I) 



Ap(;r02)- 




' for |A0| < ©s 
otherwise 



or in Fourier space 

2Ap7i(ms) 



F{£@,) 



f&s 



(56) 



(57) 



Because of the linearity of the reconstruction algorithm, the peak 
amplitude is unimportant for the shape of the radial p.s.f. We 
therefore simply set Ap - I henceforth. Unlike the discussion in 
Simon et al. (2009) for calculating the radial p.s.f., we here have 
also to factor in the tracer number density on the /th lens plane 
(and same direction 0). As this is a random variable for r{{) + 1, 
we define the p.s.f. as the radial density profile in the smoothed 
map given a matter peak 5^ (^) = F(&s{) on the /th lens plane 
marginalised over the tracer density hg (f), which is associated 
with the mass peak. This conditional mean tracer number density 
is given by 



n«(^;F)-(««(^)|5«(^) 



rJO', 



hl>F{@J)r^\e)b^\e) , (58) 

where the conditional ensemble averaae 



{x\y}n = 



J dx P(x, y)x 
J dx P(x, y) 



(59) 



is taken over all realisations of the tracer density field, and P{x, y) 
denotes the bivariate p.d.f. of the tracer number density x and 
the matter density y. The expression on the r.h.s. in the Eq. (58) 
is exact only for Gaussian statistics, which is assumed here as 
lowest-order approximation. For a differing statistics, such as a 
log-normal tracer density field (Coles & Jones 1991), we have 



to expect deviations from this expression. Evidently, the condi- 
tional average will vanish, if the correlation factor is r*''(/') = 0: 
the average tracer number density about a mass peak vanishes in 
this case. 

According to this definition, the radial p.s.f. equals the aver- 
age sight-line density profile (Eq. 77 of Simon et al. 2009): 



r- 1 r°° dee ~ ~ r. ^ 

f(S„(0,), «g(©,)] - J — FitQ,) W, [6Ui), «g(^)] (60) 



where the Wiener filter Wf is given in Eq. (55), and the vectors 
6m{e) and «g(^) {Nip elements) vanish everywhere except in their 
/th element that equals 1 and hg r'-'\f)b^'\(), respectively. The 
elements of the vector ^,n(©s) encapsulate the radial p.s.f. of the 
matter map, and the radial p.s.f. of the tracer number density map 
in the case of «g(©s). The former is the focus in the following. 

The resulting p.s.f. of the fiducial survey and a pixel size of 
©s = 1 arcmin is depicted in Fig. 2 for the cases r(0 = 0,0.8. 
Owing to the Wiener smoothing, the mass peaks are generally 
radially smeared, and their amplitude is suppressed, in particu- 
lar for very small and high redshifts. Compared to the lensing- 
only technique {r(e) = 0), however, adding tracers with r ^ 
to the map making process, clearly improves the p.s.f.: the ra- 
dial profiles are narrowed and more pronounced, the amplitudes 
are less suppressed. As the S/N in the map is improved as well, 
see left panel of Fig. 1, the redshift accuracy in the map also 
improves. The redshift accuracy is basically determined by how 
well the smoothed profiles of mass peaks at different redshifts 
can be distinguished considering the noise in the maps (see Sect. 
2.6 in Simon et al. 2012). More importantly, the positions of 
smeared peaks now coincide with the redshifts of the peaks in 
the un-smoothed map. Therefore, a moderate mix of 3D lens- 
ing and galaxy clustering information potentially lifts the z-shift 
bias. This is explored in more detail by Fig. 3, where we find 
essentially no z-shift bias for r(£) > 0.4. 
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4.4. Mock data 3D reconstruction 

As a demonstration of the mapping algorithm for a concrete 
3D matter and galaxy number density field, we apply the new 
methodology to simple mock data with survey and binning pa- 
rameters as described in Sect. 4.1. For the mocks, we place 
nine galaxy cluster scale mass peaks with a singular isothermal 
sphere (SIS) density profile and mass M200 - 6.6 x IO'^^/i^'Mq 
(SIS velocity cry = 10^ kms"' at z = 0) on a regular angular grid 
that covers Afov = 1x1 deg- with 256 x 256 pixels. The peaks 
are situated at the lens plane redshifts z = 0.15, 0.25, etc.. For 
more details, we refer the reader to Simon et al. (2009) where 
similar mock data was produced. 

We convert the model SIS mass density contrast on the lens 
planes, 6^, into pixelised shear fields 7 = Py^QS^ according to 
Sect. 2.1. Then, we randomly draw source galaxy positions 6j 
and for each source random intrinsic ellipticities e' to obtain the 
final mock catalogue source ellipticity y^'^Oj) + e'. The number 
of sources of each source redshift bin is determined by the total 
source number, fiAfov ~ 10^, and the redshift distribution p^iz). 

Every lens plane grid S^ serves as input for a constrained 
realisation of the tracer number density wj, . As the correlation 
between matter and galaxy number density is not supposed to be 
perfect, here r*''(0 - 0.8, we have to simulate some additional, 
statistically independent scatter of the tracer density about the 
matter density. For the scope of this demonstration, we assume 
Gaussian statistics for the tracer density contrast 



,.(') 



W = 



^(0 



(61) 



(0 



In this simplified case, a constrained realisation of k„ is easily 



r.(')/ 



done in Fourier space for every Fourier coefficient k), {{) inde- 
pendently. The Fourier coefficients f{t) of a Gaussian random 



Fig. 4. Example of a constrained realisation of the tracer density con- 
trast /fg on a pixelised lens plane at redshift z = 0.45. The size of the 
figure is 1 X 1 deg- with 256 x 256 pixels. At the centre is a simulated 
mass peak with SIS density profile and M200 = 6.6 x 10''* Mq/;"'. The 
adopted correlation between tracer number density and matter density 
is r(C) = 0.8. The bias factor is fo(0 = 1. 



field are independent from each other. The bivariate p.d.f. of the 
Fourier coefficients of the tracer and matter ffeld is 



P(^g,5„) 



1 



Ino-^b Vl -r2 



exp 



Kl/b^ + 'dj-lrK^Jb '' 
2cr2(l-r2) 



,(62) 



where cr^ = (5^) is the variance of the matter density ffeld, [b, r] 
are the bias parameters, and all means {Jig) - {6m) - vanish. 
In this context, by {6^, Kg} we mean either the real or imaginary 
part of the Fourier coefficients. The imaginary and real parts are 
realised independently of each other. The conditional p.d.f. of Jig 
for a given 6m, 



p{hm 



P{K„6m) 



P{Sm) 

is therefore also a Gaussian, namely with mean 

y<g\6m) = I dKgP(Kg\6m)i<g^br6m 

and variance 

cr (A-g|5in) = I di(gP{Kg\6m)i(l - b^l - r^o- . 

The variance in 6m is 

2 Ps(^ 



lA 



(63) 



(64) 



(65) 



(66) 



fov 



where Ps{i) is the power spectrum of the matter density contrast 
on the lens plane, here determined by the matter distribution of 
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Fig. 5. Impact of clipping values in the realisations of tracer number 
density fields, as demonstrated by a particular realisation of the lens 
planes redshifts z = 0. 15, 0.95. The bias and correlation factors become 
scale-dependent, most notably at low redshift, as opposed to the desired 
constant fiducial values of b(f) = 1 and r(£) = 0.8 (black horizontal 
lines). The relative error is below ~ 10% for { < 5 x 10'* = 0.4 arcmin. 



the SIS. A realisation of k„{{) from a given Sn^it) is hence ob- 
tained by 



1. multiplying djt) by biO ^Jl - r{_(f, which gives f{{); 

2. randomising the phases of the former product f{t); 

3. and adding b({)r(€)6^(e) to /(f) to yield K„{t)- 

The entire grid of kgil) obtained this way is then Fourier trans- 
formed into Kg{6), taking into account the constraint Kg{t) - 
Ki—l) for a real-valued tracer number density contrast. 

In the final step, we convert every realisation of Kg{9), such as 
the one in Fig. 4, into galaxy catalogues by drawing A^ = n„Aiov 
random galaxy positions. For thus purpose, a uniformly drawn 
position 0j is accepted, if 1 + K^iOj) > (1 + 'Cg,max)^rnd, where 
fg.max is the maximum of Kg{ff) and x^d 6 [0, 1] is an additional 
random number from a uniform p.d.f. As redshift of a galaxy, 
we choose the lens plane redshift zi and add a random Gaussian 
redshift error with r.m.s. variance o"j(z,) = 0.04(1 + zt). 

Owing to the Gaussianity of Kg in this algorithm, we also find 
pixels with 1 + Kg{9) < 0, which formally coiTesponds to neg- 
ative galaxy number densities. Since the algorithm treats such 
pixels equally as pixels with /<•;, + 1 = 0, we thereby clip all val- 
ues smaller than Kg < -\. This introduces a small bias, such that 
the tracer power spectrum, b~({)Pg{{), has not exactly the de- 
sired b{(!) = 1 on small angular scales, especially at low redshifts 
where the density fluctuations are large (Fig. 5). Furthermore, 
there is a decorrelation, r{t) < 0.8, for large £. However, the rel- 
ative error in the bias parameters is less than ~ 10% for scales 
larger than 1 arcmin, which is the angular smoothing scale in the 
following 3D maps. 

We apply the reconstruction algorithm to the mock data in 
two different ways: 

1 . In the first application, only the 3D lensing data is invoked 
for a 3D matter density map, which is equivalent to setting 
r^'\{) = for all lens planes in the Wiener filter. As tuning 
parameter, we choose a - 0.05, thus slightly more Wiener 
smoothing than in the foregoing discussion. This was nec- 
essary to attain also mass peak detections at z > 0.3 in the 
mock data. 



2. In the second run, we mix the lensing and the galaxy clus- 
tering data with known biasing parameters in order to obtain 
an updated matter density map, this time with a - 0.01 and 
j6 = 0.1. 

In both cases, the density fields on the lens planes are 
smoothed with a Gaussian kernel of ©s = 1 arcmin r.m.s. size. 
Furthermore, both maps are divided by the pixel r.m.s. noise 
variance, which we compute from 10^ noise realisations. We 
produce the noise maps by first randomising the source ellip- 
ticities and the positions of the tracer galaxies, then followed by 
rerunning the map making algorithm. The resulting S/N maps 
are shown in Fig. 6. We do not attempt to correct for the z-shift 
bias in either case. 

With lensing only (top row), we are most sensitive to struc- 
ture at low redshift, but quickly loose detection capability at 
larger radial distances, even with increased Wiener smoothing. 
This example highlights the general problems with the Wiener 
mapping, namely the radial elongation of mass peaks and the z- 
shift bias, both originating from the radial p.s.f."* In the combined 
reconstruction (bottom row), we now recover features both with 
overall higher significance and clearly reduced radial stretch, 
and also significant features at higher redshifts. As a matter of 
fact, the significance increases again at high redshifts, presum- 
ably thanks to the increasing tracer number density until z ~ 0.9. 
Moreover, the mass peaks are very close to their true redshift, 
which indicates both a reduced z-shift bias and an increase in 
the map redshift accuracy. 



4.5. Galaxy clustering noise 

The noise realisations of the data in the foregoing section are in- 
complete, as they only account for the source shape noise and the 
tracer shot-noise. In the presence of stochasticity between matter 
and tracer density, r{f) + 1 , there is a random scatter in the tracer 
density for a given matter density field, which enters the recon- 
structed matter map as additional noise component. Contrary 
to Poisson shot-noise, this noise component ("galaxy clustering 
noise" or GCN hereafter) is also present, if the number of tracer 
galaxies were infinite, and it has, in general, not a white-noise 
power spectrum. 

As can be seen in Sect. 4.4 on the generation of mock data, 
for a given matter density ~bm{t) a random component with power 
spectrum b'^(t')\\ - r^{{)]Ps{£) is added to produce a realisation 
of the tracer number density hg({), which is the GCN in the 
Gaussian case. The essential parameter for this random scatter 
is r{f), which for \r{0\ - 1 vanishes, but reaches a maximum in 
amplitude for r{£) = 0. On the other hand, a smaller r(0 also re- 
sults in a reduction of the mixing of matter and tracer density 
modes by the minimum-variance filter. On the extreme end for 
r({) - 0, the filter will not make use of any tracer information at 
all for the matter density map. 

For an assessment of the level of the GCN, we assume an 
idealised survey with Gaussian statistics and estimate the power 
of the noise that passes through the filter in comparison to the 
signal power in the mass map. The total noise power per angular 
mode in the mass map can then be approximated by 

Pn,all(0 := Pn(^) + Pgcn(0 , (67) 

* The z-shift bias in a S/N map generally differs from the z-shift bias 
in the matter density map (Simon et al. 2012). 
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Fig. 6. S/N map of the matter density field by either utilising the lensing catalogue only (top row; a = 0.05), or by mixing lensing and galaxy 
clustering data with r(f) = 0.8 and b(f) = 1 (bottom row; a = 0.01, /? = 0. 1). The right column shows a radial projection, values varying between 
0.05 and 1.15 are the lens plane redshifts. The left column panels are the corresponding projection on the sky (1x1 deg^). The circles circumscribe 
the positions of nine SIS mass peaks with M200 = 6.6 x 10'"*Mq/i"' and at increasing redshift z = 0. 15, 0.25, etc.. The shadowed areas are projected 
iso-surfaces of S/N = 3, 4, 



where Pn(^) quantifies the tracer shot-noise and source shape- 
noise inside the map, see Eq. (54), and the GCN after the appli- 
cation of the Wiener fiher Wf is 



Pgcn(^) := W,Sf "W; , 
where 



( 



diag 



0, 



,0,[«<"]McnW, 



[nf-W:^ii) 



N|n elements 



and the Gaussian GCN in the un-smoothed map is 



(68) 



(69) 



(70) 



For the previous fiducial survey and for r{f) e [0, 1 .0] and 
b{() e [0, 3.0], we find that the GCN is smaller than ~ 30% of 
the actual signal Ps in the map (Fig. 7) for angular modes larger 
than ~ 1 arcmin. For this estimate, we take the cosmic average 
Pfj {() as the signal power in the un-smoothed matter density 
field. This certainly underestimates the fluctuations in the rele- 
vant galaxy cluster regions. On the other hand, the GCN scales, 
at least for Gaussian random fields, linearly with the amplitude 
of the actual matter density fluctuations, Eq. (70), such that the 
GCN-to-signal ratio in Fig. 7 should still be a robust approxima- 
tion for Gaussian fields with more matter clustering. We find the 
GCN in the smoothed maps to be most prominent for r{C) ~ 0.6, 
declining for correlations greater or weaker than that. The GCN 
increases with the clustering b{{) of the tracers and declines for 
angular modes with wave-numbers greater than { ~ 3 x 10^ (~ 7 
arcmin). 
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Fig. 7. Estimated percentage GCN-to-signal ratio Pgcii(f)/Ps(^) Cv-axis) as function of lens plane redshift (different data points) and angular wave- 
number (x-axis). The tuning parameters are a = 0.01 and/3 = 0.1. Left: The bias factor is set to b{l) = 1, different panels correspond to different 
correlation factors r{C). Right: Varying bias factors (panels), while keeping the correlation factor at r{{) = 0.8. There is one exception in the very 
right panel that depicts a strong GCN level for r{i) = 0.2 and b{i) = 3.0. 
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every lens plane of a map, such as Fig. 6, can in principle be 
rescaled by the following recipe to incorporate the GCN. 

The pixel-to-pixel variance cr^jj in a mass map has three in- 
dependent components. 



'all 



o", + 0-- + 0-, 



shot 



(i+/gc„)f^s+t^; 



,2 
shot 



(73) 



where os is the variance in the matter density signal, o-gcn is the 
GCN variance, and cTs^ot is the source shape- and tracer shot- 
noise variance. In the former equation, we have substituted the 
GCN variance by the signal variance and the GCN-to-signal ra- 
tio /gen. A S/N map that accounts for both crsj,o, and o-gcn is 



Fig. 8. GCN variance /u^n on a pixel scale of 0^ = I arcmin radius rel- 
ative to the signal variance in the fiducial survey as function of lens 
plane redshift. The lines differ in their assumed r((). Except for one line 
(indicated), the bias factor of the tracers is b(£) = I. 



For each lens plane, we translate these /"-dependent GCN lev- 
els to the noise variance on the map pixel scale by virtue of the 
integral (Simon et al. 2009) 



[<5^cn(0s)]' 



Jo 27r 



F(^0s)r[P.gcn(O],7 



(71) 



and likewise for the signal power [Ps(/)],7 in Eq. (53) to compute 
([(Js'*(0s)]-). Fig. 8 shows the resulting GCN-to-signal ratios 
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7 gen 



\ 
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(72) 



<0 



for pixels in our fiducial mass map. Typical figures for /g(,„ 
are below 30%, but can go above this level for strongly clus- 
tered tracers. The exact value of /g^n depends on the survey and 
Wiener filter parameters, which can be determined by the above 



Ai) 



formalism. Once /g^n has been determined, the S/N levels in 
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where <5m,est(^)/o"shot on the right-hand side (rh.s.) is the S/N 
invoking shot-noise only, such as the bottom panel of Fig. 6. 
For the correction factor inside the brackets, the signal variance 
(Tshot can be estimated by employing Eq. (71) with an appro- 
priate Wiener-filtered signal power spectrum. In addition, the 
shot-noise variance, cr^, is determined by the pixel variance in 
the noise realisations, or also by Eq. (71) with the Wiener filter 
noise power spectrum Pn(0 inside the integral. For a signal vari- 
ance (Ts <& cTshot, the correction factor is roughly unity, which 
is always the case for a cosmic average matter density power 
spectrum and our fiducial survey. 

As we are mainly targeting galaxy cluster regions with lens- 
ing cartography, however, a fiducial value of as of higher vari- 
ance than in a cosmic average is likely. To obtain a more rea- 
sonable fiducial value, we construct for o-f^ an alternative signal 
power spectrum, based on the assumption of (i) randomly scat- 
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tered haloes with average number density n^is and (ii) haloes with 
a SIS matter density profile^ (Simon et al. 2009) 



€(^) = 



87r2 



3nM+z(Xi)] 



(vF 



Di 



1 



hixd^Xi ( 



(75) 



in Fourier space and SIS velocity cry. Therefore, the matter 
power spectrum for the /th lens plane is described by (Appendix 
A) 



Pf(i) = l€(^)l'«sis 



(76) 



which we insert into Eq. (53) and Eq. (71) to calculate cts- In 
accordance with our mock data, we use Wsis - 1 deg"^ and 
cr^ - lO^'kms"', which actually is a conservative upper limit 
since the SIS velocity of the mass peaks in the simulation de- 
creases with redshift due to the fixed M200 mass. We find that 
for r{{) = 0.8 and b(€) = 1 the S/N levels in Fig. 6 have to be 
reduced to ~ 75%, 90% at z = 0.15,0.25 and to values > 93% 
otherwise. These figures are roughly lower limits for r{€) G [0, 1] 
and b{() e [0, 3]. This qualitatively matches our observation in 
the mock data that there are only extra features at low redshift in 
the combined reconstruction. 



4.6. Mismatched biasing parameters 

In the foregoing discussion, we have always assumed that the 
Wiener filter devises a set of galaxy biasing parameters that ex- 
actly matches the biasing inside the data. More generally, how- 
ever, we may also use biasing parameters for the Wiener smooth- 
ing (Step-2 in Eq. 41), 



B 



:= (1 + SN„;J"'SN, 



(77) 



that are different from the biasing in the data (mismatched), as 
long as we devise for the forgoing analysis of the radial p.s.f. 
and the S/N map the correct parameter set for the actual galaxy 
biasing in the data; we would thus have different values for 
[b{{), r{€)] in Wf and Sf in Eq. (53). 

In a practical application, we may deliberately decide in 
favour of mismatched values of r{€) by choosing a value in the 
filter B just high enough to apply sufficient mixing of matter and 
tracer density modes to lift the z-shift bias. 

As extreme example for the impact of wrong bias parame- 
ters, consider the case that we adopt r{E) - 0.8 for the Wiener 
filter and assume the same value for the data, although we have 
no correlation between the matter and galaxy distribution in the 
data to which the Wiener filter applies, i.e., {KglSm) = instead 
of {kglSm) - brSm (Sect. 4.4). If this is ignored, we will arrive 
with Eq. (60) at the conclusion that the right panel of Fig. 2 is 
the p.s.f. inside the matter map instead of a p.s.f. very similar 
to the left panel (not shown), which exhibits a clear z-shift bias. 
Moreover, we will underestimate the level of GCN and hence 
overestimate the S/N levels in the map, which will lead to an 
increase of false positive detections in the map. 

4. 7. Non-Gaussianities 

For the GCN, a great uncertainty is the true nature of the galaxy 
biasing, which up to here has been treated as Gaussian, such 



^ The formalism in Sect. 2.1 is built around the mean density con- 
trast (5,^ (0) within a volume slice of width Ax, whereas the SIS density 
contrast is commonly written as integrated contrast along a l.o.s.; the 
difference between both definitions is the factor Ax- 



that a linear stochastic biasing model is a complete description 
of galaxy biasing. Although this is probably a fair description 
for large smoothing scales beyond ~ 8/i"'Mpc, it is increas- 
ingly inaccurate at smaller scales where the biasing is expected 
to become non-linear (e.g., Yoshikawa et al. 2001). As imme- 
diate consequence, the Gaussian GCN analysis above probably 
underestimates the clustering noise on small angular scales. The 
detection rate of false positives in the mass map would rise as a 
consequence. 

Moreover, non-linearities in the biasing may have unwanted 
side-effects, such as a mass map being strongly skewed towards 
the galaxy distribution. This may be illustrated by the high-peak 
bias model (Kaiser 1984; Bardeen et al. 1986), which is non- 
linear and deterministic. In this model, galaxies are found inside 
peaks of the matter density field that exceed a certain density 
threshold ("high-peaks"); outside those peaks, the galaxy num- 
ber density is exactly zero. This model is deterministic, as it as- 
sumes the same Kg for a given S^- Nevertheless, the correlation 
factor is not exactly r({) - 1, since the galaxy number density 
always vanishes outside the high-peaks, even though the matter 
density may be varying. Due to r{() + 0, there is a linear mixing 
of tracer and matter density modes in the Wiener filter which, 
however, suppresses the matter density outside the high-peaks, 
since there we have k^ - -1. The mass map outside the high- 
peaks is probably underrepresented as a result. We have shown 
that the Wiener filter is biased by construction, as it is applying 
smoothing to the maps. But this example shows that the exact 
appearance of the bias also depends on the details of the galaxy 
biasing mechanism. It is likely that for a more accurate mapping 
algorithm on non-linear scales, a more elaborate modelling of 
the statistics and the detailed galaxy biasing is required, perhaps 
in the vein of Jasche et al. (2010) but also incorporating lensing. 

5. Summary and conclusions 

In this study, we introduce and discuss a new technique for the 
synergy of lensing and galaxy clustering information in order to 
map the 3D matter density. This technique devises a minimum- 
variance filter based upon both the 3D shear information from 
radially binned sources and the angular distribution of a coarsely 
binned galaxy tracer population with possibly inaccurate red- 
shifts. The filter essentially proceeds in two steps: Step-1 pro- 
duces an unbiased estimate of the matter and tracer number den- 
sity, while Step-2 linearly combines both to produce smoothed 
density fields with an enhanced S/N. The estimator becomes bi- 
ased through Step-2, but it is consistent as it becomes asymptoti- 
cally unbiased when the S/N in the data approaches infinity. The 
smoothing of the maps is uniquely defined in Eq. (77), which 
and can be applied to theoretical maps of the matter and galaxy 
number density for a quantitative comparison to the data. 

As novelty, the filter linearly mixes the matter and galaxy 
(tracer) number density with a weighting that depends on the 
galaxy bias of the tracer population. In principle, the biasing pa- 
rameters can be constrained from the data itself by using lensing 
techniques (cf., Hoekstra et al. 2002; Simon 2012). In a situa- 
tion, where the tracers have no statistical correlation to the mat- 
ter density field, there is no mixing and hence no gain from a 
synergy. The level of smoothing and mixing is regulated by two 
independent tuning parameters that rescale the noise covariance 
inside the filter. 

For the discussion of the filter, we assume a fiducial survey 
with mean redshift I = 0.84, a source number density of « = 30 
galaxies per arcmin^, in accordance with Simon et al. (2009), 
and a tracer number density of «g = 3 galaxies per arcmin^ that 
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is almost unbiased but has some stochasticity (Sect. 4.4). Most of 
the information on the matter density originates from the clus- 
tering of the tracer population, at least when the biasing is ex- 
actly known. We have chosen to keep the influence of the tracer 
density maps on the mass map at a low level by artificially re- 
ducing the S/N level of the tracers in the Wiener filter, which is 
done by appropriately choosing the tuning parameters (a - 0.01, 
/? = 0.1). Moreover, in future applications, the mixing could ad- 
ditionally be reduced by adopting a correlation factor inside the 
Wiener filter that is smaller than the actual correlation in the data. 

The S/N in the matter density map is improved by the syn- 
ergy, especially at higher redshifts. However, at our adopted level 
of mixing, significant detections are basically still restricted to 
peaks of a galaxy cluster mass scale. By quantifying the aver- 
age radial smoothing kernel orp.s.f. in the map, we have shown 
that the most useful improvement of a synergy is a removal of 
the z-shift bias, which is known from lensing-only methodolo- 
gies based on linear filters. Even for relatively moderate levels 
of correlation between the matter and tracer density, r > 0.4, the 
z-shift bias is removed from the mass map. Likewise, the peaks 
in the mass map are also less radially smeared out owing to the 
narrower radial p.s.f., in particular for high correlation factors. 
An increase of the S/N and less smearing at the same time also 
results in a higher redshift accuracy of the mass peaks. 

The synergy introduces a new kind of noise into the mass 
maps (galaxy clustering noise; GCN) that stems from a random 
scatter of the galaxy number density about a fixed matter den- 
sity field (stochastic biasing). We derive estimates for the GCN 
in the Gaussian regime and show that it can relatively easy ac- 
counted for through rescaling the S/N maps that are obtained by 
randomising the source shapes and the galaxy tracer positions. 
For the fiducial survey discussed in the paper, this correction 
mainly aff'ects the lens planes with z < 0.3. 

The GCN discussion, however, assumes Gaussian statis- 
tics which only applies on large smoothing scales greater than 
~ 8 /i^'Mpc. On smaller scales, the statistics is likely to be more 
complex, possibly resulting in spurious false detections in the 
mass map whenever a Gaussian GCN model is used. Moreover, 
we argue that the detailed bias of the filter depends on the de- 
tails of the galaxy biasing mechanism in the non-linear regime. 
More elaborate models or computer simulations are needed to 
investigate this eflect. 
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Appendix A: SIS power spectrum 

A set of A^h haloes with positions 6j on the /th lens plane, and 
an average, axial-symmetric matter density contrast ^sisCI^I) pro- 
duces the combined density contrast 
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where the second equation on the r.h.s. is the Fourier transform 
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of 6,^(6). Averaging the two-point correlator of the density in 
Fourier space over all halo positions, results in 
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We ignore the clustering of the haloes over the field-of- view Afov , 
so that the two-halo term in the second sum vanishes, and 
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= (27r)25D(/'i+^2)Ss.s(A)5sis(^2)nsis (A.4) 

^ {2nf6^{{i+{2)Pf((x), (A.5) 

where ;isis '■- A^hMfov expresses the mean number density of 
haloes. Therefore, we obtain in this scenario 
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